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Abstract. Using the axially-symmetric time-dependent Gross-Pitaevskii equation we study the phase co- 
herence in a repulsive Bose-Einstein condensate (BEC) trapped by a harmonic and an one-dimensional 
optical lattice potential to describe the experiment by Cataliotti et al. on atomic Josephson oscillation 
[Science 293, 843 (2001)]. The phase coherence is maintained after the BEC is set into oscillation by a 
small displacement of the magnetic trap along the optical lattice. The phase coherence in the presence 
of oscillating neutral current across an array of Josephson junctions manifests in an interference pattern 
formed upon free expansion of the BEC. The numerical response of the system to a large displacement of 
the magnetic trap is a classical transition from a coherent superfluid to an insulator regime and a subse- 
quent destruction of the interference pattern in agreement with the more recent experiment by Cataliotti 
et al. [e-print cond-mat/0207139 . 

PACS. 03.75.-b Matter waves - 03.75.Lm Tunneling, Josephson effect, Bose-Einstein condensates in peri- 
odic potentials, solitons, vortices and topological excitations - 03.75Kk Dynamic properties of condensates; 
collective and hydrodynamic excitations, superfluid flow 



1 Introduction 

The observation of an oscillating Josephson current across 
the boundaries of a one-dimensional array of potential 
wells, usually generated by a standing-wave laser field 
and commonly known as an optical lattice potential, in a 
trapped cigar-shaped Bose-Einstein condensate (BEC) by 
Cataliotti et al. is a clear manifestation of macroscopic 
quantum phase coherence. So far the Josephson effect has 
been confirmed in superconductors with charged electrons 
and in liquid helium [3]. 

The recent experimental observation of BEC in trapped 
alkali-metal atoms has offered new possibility of the 
confirmation of Josephson effect in neutral quantum flu- 
ids with an array of quasi one-dimensional Josephson junc- 
tions not realizable in superconductors. The experimental 
loading of a cigar-shaped BEC in both one- and 
three-dimensional [J| optical lattice potentials has allowed 
the study of quantum phase effects on a macroscopic scale 
such as interference of matter waves jBJ- There have been 
several theoretical studies on different aspects of a BEC in 
a one- [5] as well as three-dimensional JU] optical lattice 
potentials. The phase coherence between different sites of 
a trapped BEC on an optical lattice has been established 
in recent experiments 00H51IZI through the formation of 
distinct interference pattern when the traps are removed. 

Cataliotti et al. [Till 1| have provided a quantitative 
measurement of the formation and evolution of interfer- 
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ence pattern upon free expansion of a cigar-shaped trapped 
BEC of repulsive Rb atoms on an optical lattice and har- 
monic potentials after the removal of the combined traps. 
The phase coherence in a BEC trapped in a standing- wave 
optical-lattice is responsible for the formation of a distinct 
interference pattern upon free expansion as observed in 
several recent experiments 1 1 •") 7 1 1 . Cataliotti et al. 
PQ also continued their investigation to a BEC oscillating 
on the optical lattice via quantum tunneling and found 
that the phase coherence between different sites is main- 
tained during oscillation initiated by a sudden shift of the 
magnetic trap along the optical axis. 

The phase-coherent BEC on the optical lattice is a 
quantum superfluid |7] and the atoms in it move freely 
from one optical site to another by quantum tunneling. 
However, the classical movement is prohibited through the 
high optical potential traps. It has been demonstrated in 
a recent experiment by Greiner et al. [J] that, as the op- 
tical potential traps are made much too higher, the quan- 
tum tunneling of atoms from one optical site to another 
is stopped resulting in a loss of superfluidity and phase 
coherence in the BEC. Consequently, no interference pat- 
tern is formed upon free expansion of such a BEC which 
is termed a Mott insulator state. This phenomenon repre- 
sents a quantum phase transition (with energy nonconser- 
vation in the tunneling process) and cannot be accounted 
for in a classical mean-field model based on the Gross- 
Pitaevskii (GP) equation [12] . 
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Following a suggestion by Smerzi et al. j!3| . more re- 
cently Cataliotti et al. ^3] have demonstrated in a novel 
experiment the loss of phase coherence and superfluid- 
ity in a BEC trapped in a optical-lattice and harmonic 
potentials when the center of the harmonic potential is 
suddenly displaced along the optical lattice through a dis- 
tance larger than a critical value. Then a modulational 
instability takes place in the BEC and it cannot reorga- 
nize itself quickly enough and the phase coherence and 
superfluidity of the BEC are destroyed. The resulting mo- 
tion of the condensate is not oscillatory in nature. The loss 
of phase coherence is manifested in the destruction of the 
interference pattern upon free expansion. However, for dis- 
placements smaller than the critical distance the BEC can 
reorganize itself and the phase coherence and superfluidity 
are maintained [Dll4| . Recently, a new classical mechanism 
for the loss of superfluidity of a BEC arising from a non- 
linear modulation of the scattering length has been sug- 
gested rjHl- Distinct from the quantum phase transition 
observed by Greiner et al. |7|, these modulational insta- 
bilities responsible for the destruction of phase coherence 
are classical in nature and can be described |13U14| by the 
mean-field model. Hence in the present paper we present 
a mean-field description of the experiments by Cataliotti 
et al. PHT4] to see how well and how far it can describe the 
observed phenomena. Specifically, we consider the numer- 
ical solution of the axially-symmetric CP equation rj2] m 
an optical and a harmonic trap. 

Cataliotti et al. [PTEl] provided a theoretical account of 
their study using the tight-binding approximation for the 
full wave function in the presence of the periodic optical 
potential wells. Also, there has been a preliminary attempt 
to explain some features of this experiment using one- 
dimensional mean- field models |11U16| . In reference [T?)j 
a dynamical solution of one-dimensional CP equation was 
used; whereas in reference an one-dimensional model 
of interference was developed using superposition of ana- 
lytical matter waves, which is reasonable in the absence of 
nonlinear atomic interaction. Although, the tight-binding 
approximation and these one-dimensional models could be 
reasonable for the study of some aspects of the experiment 
of Cataliotti et al. [Dll4| , here we compare the results with 
the complete solution of the three-dimensional mean-field 
Hamiltonian via the nonlinear GP equation |12| . 

In section 2 we present the mean-field model based on 
the axially-symmetric timc-dcpcndcnt nonlinear GP equa- 
tion. In section 3 we present the numerical results and 
finally, in section 4 we present the conclusions. 



2 Mean-field Model 

The time-dependent BEC wave function ^(r; t) at posi- 
tion r and time t is described by the following mean-field 
nonlinear GP equation |12j 



-if l 



dt 



2m 



+ V{v)+gN\^{v;tt 



#(r;t)=0,(l) 



where m is the mass and N the number of atoms in the 
condensate, g = 4irh 2 a/m the strength of interatomic in- 



teraction, with a the atomic scattering length. In the pres- 
ence of the combined axially-symmetric and optical lattice 
traps V(r) — ^mw 2 (p 2 + v 2 y 2 ) + V p^ where ui is the angu- 
lar frequency of the harmonic trap in the radial direction 
p, vui that in the axial direction y, with v the aspect ratio, 
and V^pt is the optical lattice potential introduced later. 
The normalization condition is J c?r|if r (r; t)\ 2 = 1. 

In the axially-symmetric configuration, the wave func- 
tion can be written as ^(r,i) = ip(p,y,t), where < 
p < oo is the radial variable and — oo < y < oo is the 
axial variable. Now transforming to dimcnsionlcss vari- 
ables p = y/2p/l, y — \/2y/l, r = tui, I = ^/h/(mui), and 



ip(p,y;T) = pJ I 3 / y/8ip(p, y; t), equation @J becomes ^3 
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p 



<P(P,V;t) = 0,(2) 



where n = Na/l. In terms of the one-dimensional proba- 
bility P(y,t) = 2ir J Q dp\ip(p,y,T)\ 2 / p, the normalization 
of the wave function is given by dyP(y,t) = 1. The 
probability P(y, t) is useful in the study of the present 
problem under the action of the optical lattice potential, 
specially in the investigation of the formation and evolu- 
tion of the interference pattern after the removal of the 
trapping potentials. 

In the experiment of Cataliotti et al. pQ with repulsive 
87 Rb atoms in the hyperfine state F = l,mp = —1, the 
axial and radial trap frequencies were vui — 2tt x 9 Hz 
and u> = 2ir x 92 Hz, respectively, with v = 9/92. The op- 
tical potential created with the standing- wave laser field of 
wavelength A = 795 nm is given by V -p^ = VqEr cos 2 (fcz,z), 

with E R = h 2 k 2 L /(2m), k L = 2tt/A, and V (< 12) the 
strength. For the mass m = 1.441 x 10~ 25 kg of 87 Rb 
the harmonic oscillator length I = y/h/ (raw) = 1.126 /im 
and the present dimcnsionlcss length unit corresponds to 
l/y/2 = 0.796 fim. The present dimensionless time unit 
corresponds to ui^ 1 = 1/(2tt x 92) s = 1.73 ms. Although 
we perform the calculation in dimcnsionlcss units using 
equation (J2J, we present the results in actual physcial 
units using these conversion factors consistent with the 
experiment by Cataliotti et al. . In terms of the dimen- 
sionless laser wave length Ao = \[2\jl ~ 1, the dimension- 
less standing- wave energy parameter Er/Qiui) = 4tt 2 /Xq. 
Hence in dimensionless unit V^p^ of equation @ is 
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We solve equation (TJJl numerically using a split-step 
time-iteration method with the Crank-Nicholson discretiza- 
tion scheme described recently ^HJ- The time iteration 
is started with the known harmonic oscillator solution 



of equation with 



0: <p(p,jj) = K(8tt 3 )] 



3 YlV4 n 



e -(p +»y )/4 with chemical potential p = (l+u/2) [T7|. For 
a typical cigar-shaped condensate with ^~0.irjQ/2~lis 
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Fig. 1. P(y,t) vs. y and t for the ground-state BEC with n = 
10 and Vo = 6 for < t < 100 ms and (a) -30 /im < y < 30 
/tm and (b) —6 /xm < y < 6 /jm. The optical potential for 
—6 (jm < y < 6 /jm is shown in (c). The harmonic oscillator 
potential is negligible on this scale. 

much smaller than the typical depth of the optical poten- 
tial wells E R /(hbj) = 47r 2 /Ag ~ 40 so that p, « E r /(Twj) 
and the passage of condensate atoms from one well to 
other can only proceed through quantum tunneling. The 
nonlinearity as well as the optical lattice potential param- 
eter Vo are slowly increased by equal amounts in lOOOOn 
steps of time iteration until the desired value of nonlinear- 
ity and optical lattice potentials are attained. Then, with- 
out changing any parameter, the solution so obtained is 
iterated 50 000 times so that a stable solution is obtained 
independent of the initial input and time and space steps. 
The solution then corresponds to the bound BEC under 
the joint action of the harmonic and optical traps. 

3 Numerical Results 

First we consider the BEC formed on the optical lattice 
for a specific nonlinearity. In the present study we take 
nonlinearity n = 10 and optical lattice strength Vo = 6 
except in figure 5 where we use a variable Vq. We con- 
sider the ground-state wave function in the combined har- 




Fig. 2. P(y, t) vs. y and t for the BEC of figure 1 after the 
removal of combined traps at t = for a lattice defined by 
p < 25 /im and (a) —320 fim < y < 320 fim and (b) —160 
/im < y < 160 /im. 



monic and optical lattice potentials. The one-dimensional 
pattern in the axial y direction is most easily illustrated 
from a consideration of the probability P(y, t) at different 
times. In figure 1 (a) we plot the frontal view of P(y,t) 
for < t < 100 ms and — 30/im < y < 30/mi. In this in- 
terval of y, there are 150 wells of the optical potential and 
as many maxima and minima in P(y, t), which cannot be 
visualized clearly in figure 1 (a). In the actual experiment 
200 wells were typically populated, which corresponds to a 
larger condensate than considered in this numerical simu- 
lation. For the limitation in computer processing time we 
had to stick to a smaller condensate. In figure 1 (b) we 
show a close-up of figure 1 (a) for —6 /im < y < 6 /im 
containing 30 wells. The corresponding optical potential 
is shown in figure 1 (c), which clearly shows the 30 wells. 
From figure 1 (b) one can count 30 maxima and 30 minima 
in probability P(y,t). 

As the present calculation is performed with the full 
wave function without approximation, phase coherence 
among different wells of the optical lattice is automati- 
cally guaranteed. As a result when the condensate is re- 
leased from the combined trap, a matter-wave interference 
pattern is formed in few milliseconds. The atom cloud re- 
leased from one lattice site expand, and overlap and inter- 
fere with atom clouds from neighboring sites to form the 
robust interference pattern due to phase coherence. No in- 
terference pattern can be formed without phase coherence. 
The pattern consists of a central peak at y = and two 
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Fig. 3. P(y,t) vs. y and t for an oscillating BEC on optical 
lattice after a displacement of the magnetic trap through 25 /jm 
along the optical axis and upon the removal of the combined 
traps at t = 35 ms (hold time). 



symmetrically spaced peaks, each containing about 10% 
of total number of atoms, moving apart from the central 
peak P"ITTT| 

The simulation of the formation of the interference pat- 
tern is performed by loading the preformed condensate 
of figure 1 on two lattices with p < 25 pro. and (a) 320 
/im > y > —320 /im, and (b) 160 pm > y > —160 pm. 
which will permit the study of the evolution of the inter- 
ference pattern on a large interval of space and time. The 
interference pattern is formed by suddenly removing the 
combined traps at time t — 0. The time evolution of the 
system is best illustrated via P(y,t) and we plot in fig- 
ures 2 (a) and (b) P(y, t) vs. y and t for lattices (a) and 
(b), respectively. The dynamics is symmetric about y = 
and only P(y,t) for positive y is plotted in figure 2 (a). 
In these plots we can clearly see the central condensate 
and the moving interference peak(s). The peaks spread 
unobservably slowly as they propagate, even after reflec- 
tion from the boundary or after crossing each other. The 
phase coherence between the components of BEC at differ- 
ent sites of optical lattice is responsible for the generation 
of the interference pattern with very little or practically 
no spreading. Without the initial phase coherence over a 
large number of lattice sites, a repulsive condensate in the 
absence of a trap will disappear in few milliseconds |19| . 
Each of the moving interference peaks is similar to atom 
laser |Hl20| which can be used in the scattering of two 
coherent BECs and other purposes. 

We have also examined the wave function ip(p, y, t) at 
different times (not reported here). We find that there is 
virtually no spreading of the wave function in the axial 
y direction during few hundred milliseconds. The phase 
coherence in the axial direction due to the optical lattice 
is responsible for the localization of the peaks. 

Next we consider an oscillating BEC in the combined 
harmonic and optical traps. If we suddenly displace the 
magnetic trap along the lattice axis by a small distance 
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Fig. 4. Center of mass positions of the three interference peaks 
of the expanded condensate after 20 ms of free expansion vs. 
hold time of the oscillating BEC. The magnetic trap is dis- 
placed through a distance of 25 /jm. The + symbols denote 
the results of simulation which are joined by full lines to show 
the correlated oscillation of the three peaks. 



after the formation of the BEC in the combined trap, the 
condensate will be out of equilibrium and start to oscil- 
late. As the height of the potential-well barriers on the 
optical lattice is much larger than the energy of the sys- 
tem, the atoms in the condensate will move by tunneling 
through the potential barriers. This fluctuating transfer of 
Rb atoms across the potential barriers is due to Joseph- 
son effect in a neutral quantum liquid. The experiment 
of Cataliotti et al. [TUTT] demonstrates that the phase co- 
herence between different wells of the condensate is main- 
tained during this mass transfer and a matter-wave in- 
terference pattern with three peaks is formed after the 
removal of the joint trap. The peaks of the expanded con- 
densate oscillate in phase, thus showing that the quantum 
mechanical phase coherence is maintained over the entire 
condensate. They studied this problem experimentally in 
some detail by varying the time of oscillation of the BEC 
(hold time) before removing the combined trap . 

To see if the present simulation can represent the es- 
sential features of the phase coherence of the oscillating 
BEC, we load the GP equation with the BEC of the com- 
bined harmonic oscillator and optical traps on a lattice 
defined by 200 /im > y > —200 /im and p < 25 psm and 
suddenly displace the harmonic trap along the optical axis 
by 25 /im. The BEC starts to oscillate and we allow the os- 
cillation to evolve through a certain interval of time, called 
hold time, before the removal of the combined traps. The 
interference pattern is observed after some time of free 
expansion and the positions of the interference peaks are 
noted. In figure 3 we plot the one-dimensional probability 
P(y, t) vs. y and t after an initial evolution of the oscilla- 
tion during 35 ms and observe the interference pattern for 
160 ms. In this case, unlike in figures 2, the large central 
peak does not stay at rest and the sizes and positions of 
the two smaller peaks are not symmetrical around y = 0. 
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Fig. 5. The frequency of the atomic current in the array of 
Josephson junctions as a function of optical lattice strength: • 
with error bar — experiment of Cataliotti et al. pQ; A — tight 
binding calculation [Q; * — present calculation. 



The clear formation of the interference pattern with very 
little spreading even after reflection at the boundaries and 
its propagation for more than 160 ms is noted in figure 3 
which confirms the phase coherence. 

To study the phase coherence in detail we plot in fig- 
ure 4 the positions of the expanded interference peaks 
after 20 ms of free expansion for different hold times of 
oscillation in the displaced harmonic potential. We find 
that the interference peaks oscillate in phase showing the 
phase coherence in the oscillating BEC. Similar oscilla- 
tion was also observed in the experiment of Cataliotti et 
al. . From figure 4 we find that the period of this oscil- 
lation is about 170 ms corresponding to a frequency of 5.9 
Hz, which is very close to the experimental result exhib- 
ited in figure 3 of reference [Q . To make a more complete 
comparison with figure 3 of reference £Q we calculated 
the frequency of atomic current in the array of Joseph- 
son junctions for different Vq and the results are shown in 
figure 5 where we plot the present frequencies as well as 
those of the experiment of Cataliotti et al. £Q and of their 
tight-binding calculation. From figure 5 we see that the 
complete solution of the GP equation has led to results 
in agreement with the experiment of Cataliotti et al.. The 
agreement of the present calculation in figure 5 performed 
with a smaller condensate with experiment demonstrates 
that the frequency of atomic current is mostly determined 
by the strength of the optical lattice strength and is rea- 
sonably independent of the size of the condensate. 

Finally, we consider the destruction of superfluidity 
in the condensate when the center of the magnetic trap 
is displaced along the optical lattice by a distance larger 
than the critical distance and the BEC is allowed to stay 
in this displaced trap for an interval of time (hold time). In 
this case the BEC does not execute an oscillatory motion 
but its center moves very slowly towards the new center 
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Fig. 6. P(y,t) vs. y and t for a BEC on optical lattice after 
a sudden displacement of the magnetic trap through 120 (im 
along the optical axis and upon the removal of the combined 
traps after hold times (a) 0, (b) 35 ms, and (c) 70 ms. The 
time evolution is stopped upon 27.8 ms of free expansion after 
the removal of the combined traps. 



of the magnetic trap. The destruction of superfluidity and 
phase coherence for a larger hold time in the displaced 
trap manifests in the disappearance of the interference 
pattern upon free expansion as noted in experiment |14|. 
As in that experiment, we consider a displacement of the 
magnetic trap through 120 /im and allow the condensate 
to freely expand for 27.8 ms after different hold times in 
the optical and displaced magnetic traps. 

For numerical simulation we load the BEC of figure 1 
on a lattice with p < 25 /im and 200 (im > y > —200 /im 
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and study the its evolution after an initial displacement 
of 120 /im of the magnetic trap for hold times 0, 35 ms, 
and 70 ms. As in the experiment no oscillatory motion of 
the BEC is noted in the displaced trap. The corresponding 
probability densities are plotted in figures 6 (a), (b), and 
(c), respectively. For hold times and 35 ms prominent 
interference pattern is formed upon free expansion as we 
can see in figures 6 (a) and (b). In these cases three sepa- 
rate pieces of interference patterns corresponding to three 
distinct trails can be identified. However, as the hold time 
in the displaced trap increases the maxima of the inter- 
ference pattern mixes up and finally for a hold time of 70 
ms the interference pattern is completely destroyed as we 
find in figure 6 (c) in agreement with the experiment |14| . 

As the BEC is allowed to evolve for a substantial inter- 
val of time after a large displacement of the magnetic trap 
along the optical axis a dynamical instability of classical 
nature sets in and the system can not evolve maintain- 
ing the phase coherence |13U14j . This has been explicitly 
demonstrated in the present simulation which results in 
the destruction of the interference pattern. 



4 Conclusion 



In conclusion, to understand theoretically the experiments 
by Cataliotti et al. [TltH] . we have studied in detail the 
phase coherence along a cigar-shaped condensate loaded 
in a combined axially-symmetric harmonic trap and op- 
tical lattice trap using the solution of the mean-field GP 
equation. Upon removal of the combined traps, the forma- 
tion of an interference pattern clearly demonstrates the 
phase coherence over a very large number of optical lat- 
tice sites. Each of the moving interference peaks formed of 
coherent matter wave is similar to a atom laser observed 
experimentally (4"ll20| . The phase coherence along the op- 
tical lattice axis of the condensate is maintained even if 
the initial BEC is set into oscillation by suddenly shifting 
the harmonic trap along the optical axis through a small 
distance and keeping the BEC in the displaced trap for 
a certain hold time. This is clearly demonstrated by the 
noted correlated oscillation of the condensate peaks after 
free expansion for different hold times. The present mean- 
field model provides a proper account of the frequency 
of atomic current in the array of Josephson Junctions in 
agreement with experiment £Q. 

However, if the initial displacement of the harmonic 
trap along the optical axis is larger than a critical value 
and the BEC is maintained in the displaced trap for a cer- 
tain time, the phase coherence is destroyed. Consequently, 
after release from the combined trap no interference pat- 
tern is formed in agreement with experiment |14|. 



The work was supported in part by the CNPq and FAPESP 
of Brazil. 
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